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Dynamical Susceptibility in KDP-type Crysals above and below T c II 
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The path probability method (PPM) in the tetrahedron-cactus approximation is applied 
to the Slater- Takagi model with dipole-dipole interaction for KH^PC^-type hydrogen-bonded 
ferroelectric crystals in order to derive a small dip structure in the real part of dynamical 
susceptibility observed at the transition temperature T c . The dip structure can be ascribed to 
finite relaxation times of electric dipole moments responsible for the first order transition with 
contrast to the critical slowing down in the second order transition. The light scattering inten- 
sity which is related to the imaginary part of dynamical susceptibility is also calculated above 
and below the transition temperature and the obtained central peak structure is consistent 
with polarization fluctuation modes in Raman scattering experiments. 

KEYWORDS: KDP(KH 2 P0 4 ), Slater-Takagi model, PPM(path probability method), dynamical susceptibility, light 
scattering intensity 
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1. Introduction 

In the previous paper 1 ) (this article is referred to as 
I henceforth), we have succeeded in calculating the 
dynamical susceptibility for the Slater-Takagi model 2 ) 
above and below the transition temperature T c by mak- 
ing use of the analytical solution for spontaneous polar- 
ization in the tetrahedron-cactus approximation of the 
cluster variation method(CVM). 3 ) This approximation 
is the basic one in the CVM for taking into account the 
ice-rule for protons of KDP-type hydrogen bonded ferro- 
electrics. However, since the Slater-Takagi model always 
yields the second order phase transition, the dynamical 
susceptibility at T — T c almost vanishes as a result of 
critical slowing down of the most dominant mode gov- 
erning the transition. This result is different from the 
experimental data of KDP (KH2PO4) which show only 
a small dip at T — T c in the real part of dynamical sus- 
ceptibility versus temperature graph. 4 ) 

It is believed that KDP undergoes the first order tran- 
sition close to the second order. Senko 5 ) phenomeno- 
logically introduced the dipole-dipole interaction to the 
Slater-Takagi model and Silsbee et a/. 6 ' revealed that the 
phase transition of this model with dipole-dipole inter- 
action changes from the second order to the first order 
transition according to increase of the strength of dipole- 
dipole interaction. Based on this model Matsumoto and 
Ogura 7 ) have calculated the light scattering intensity re- 
lated to the imaginary part of the dynamical susceptibil- 
ity in the tetrahedron-cactus approximation. However, 
their calculation seems to be limited to the paraelectric 
phase. In the present manuscript we would like to show 
the small dip structure at the transition temperature in 
the experimental data of the real part of dynamical sus- 
ceptibility versus temperature for KDP-type crystal. We 
will also calculate the light scattering intensity below 
the transition temperature as well as above it. For that 
purpose we apply the tetrahedron cactus approximation 
of the PPM (Path Probability Method) 8) to the Slater- 
Takagi model with phenomcnologically induced dipolc- 
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dipole interaction since we expect that the most domi- 
nant mode governing the transition brings about a finite 
relaxation time at the transition temperature. 

2. Formulation 

There are iV PO4 tetrahedra and 27V protons in KDP- 
type crystals as shown in Fig. 1. The pseudo-spin Ising 
Hamiltonian Ti for a configuration of 27V protons has a 
form 



'H = ^ [ Wo(Cj, (Tj, CTfe, (T/) 

(ijki) ^ 



fa + (Tj + (Tfe + <7l) 



(1) 



with 



and 



Ho((Ti,aj,ak,(Ti) 

- V 4 a i a j cr k cri - V 5 (cr i a k + UjQi) + C, (2) 

Hdipoio = - 2J Vd 2 m (<7j + (Tj + a k + <Ji) 

(ijkl) (3) 

+ 27VA^ d 2 m 2 , 

where the sum (ijkl) runs over four protons i,j,k,l 
around each PO4 tetrahedron in the crystal, fj,^ is the 
magnitude of an electric dipole moment associated with 
a complex K-H2PO4, E is an external electric field, A is 
a parameter representing long range dipolar sum in the 
phenomenologically introduced dipole-dipole interaction 
^dipole and m is the electric polarization defined by ther- 
mal average over Ti as 



(0i) = (oj> = (o-fc) = (07) 



(4) 



and 27VA/id 2 TO 2 in 7idipoio is taken care of over-counting 
due to molecular field. 
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Fig. 1. The projection of atomic arrangement of KDP-type crys- 
tals on (001) plane. The variables <r; , <Tj , a k , <t; show the four 
different pseudo-spins for protons around a PO4 tetrahedron. 
The number described in a PO4 tetrahedron represents a rela- 
tive height of a PO4 tetrahedron. 
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Fig. 2. Configuration energy (£2 > £1 > £0 > 0), alloted dipole 
moment along c-axis and occurrence probability of proton con- 
figuration around PO4 tetrahedron. 



As is seen in Fig. 1, we use a convention that when the 
i-th proton is located on the closer site to an O atom at 
the bottom (top) of the PO4 tetrahedron along the easy 
c-axis, the i-th proton takes <Ji = +1(— 1). The energy 
parameters V 2l V4, V5 and C are connected to the energy 
parameters in the Slater- Takagi model shown in Fig. 2 
as 



V 2 =e 2 /8, Vi = -£ /4 + £i/2 -e 2 /8, 

Vb = £o/4 - £ 2 /8 , C = £o/4 + £l /2 + £ 2 /c 



(5) 



As shown in I, we now apply the path probability 
method (PPM) 8 ) in the tetrahedron-cactus approxima- 
tion to the present system in order to find the time evo- 
lution equation of the system with various proton con- 
figurations. After some manipulations of the PPM we 
can make a generating function from which a set of time 
evolution equations is derived through differentiation of 



an interaction parameter set L = {L\, L 2 , £3, L4, L5). 
The generating function is given by 

g(L) = -Trp 1 (<j l ,t)e- 2L ^ 



Tr 

jkl 



Pi(o"»,t) 



AiHo(cri,aj,ak, cri) 

= Ho(—cri,(jj,ak,(7i) — Ho((Ti,aj,crk,(7i), 

Li = -(l3^ d E + 2(3\^ d m(t)) , 

L 2 =(3V 2 , L 3 =pV 3 , U = pV 4 , L 5 =(3V 5 



(6) 



where Tr^ and Trjki denotes a trace operation _ ±1 
and £^,^,^=±1' respectively, (3 = l/fc B T is the 
inverse temperature with Boltzmann's constant fce, 
AiHo(<Ji,<Jj,&k,vi) defines an energy change under 
transfer of i-th proton between the double well potential 
on a hydrogen bond and To is its microscopic relaxation 
time of an isolated proton. When A = 0, the generating 
function Q(L) coincides with that of I. Thus, the effect 
of dipole-dipolc interaction modifies only the external 
electric field E by E + 2A[i < im(t). Further, state vari- 
ables pi(ai,t) and Pi(ai, aj, a k , 07, t) represent, respec- 
tively, the probability of finding the site Cj of a proton 
in the i-th hydrogen bond at time t and the probability 
of finding the sites <7j , aj ,(Jh,cri of protons i,j,k,l around 
a PO4 tetrahedron at time t. The state variables pi and 
P4 are given by 

Pi{a il aj,a k ,ai,t) 
24" 



1 / mi(t) , . 



m 2 {t) 

H ^ ( <7 i (T j + a j a k + Ok^l + CT;(7i) 

m 3 (t) 

H ^ {PiPjPk + Vj&k&l + O k (TlOi + (7l<Ti(Jj) 

+ mi{t)(Ji(jj(Jk(Ti H — \<J l a k + <jjai) 

(7) 

Here let us introduce the Ising spin space vector er 
around a PO4 tetrahedron as 



/ 



&i + (Tj + (Jk + &l 
(TiCFj + (JjCJk + &k&l + CTjCTi 
<Ji<Jj<7 k + CFjCF k ai + O k OlOi + Gl<JiGj 



\ 



V 



(Ti<Tj<J k <Jl 
O l O k + <Jj<Jl 



(8) 



/ 



Then, in the tetrahedron-cactus approximation of the 
PPM, the homogeneous state of the present system at 
time t is found to be described by the five order param- 
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eters defined by 
( 



(mi\ 

m 2 
mz 
m 4 
\m 5 J 



(<Ti + (Tj + CTfe + ai) \ 
(<7i(Tj + (JjCFk + CTkCTl + OlOi) 

(aitjjtJk + Ujdkcri + (TkcriUi + oiOiOj) 
(aiCTjCrkai} 



(9) 



where each order parameter represents the correlation 
of protons i, j, k, I around a PO4 tetrahedron at time t, 
(• • • ) t is a thermal average at time t and m(t) = mi(i)/4 
is the long range order parameter due to electric polar- 
ization per proton. Then a set of kinetic equations for 
hve order parameters is given in a convenient form as 



dmj(t) 
dt 



4 lim 



dQ{L) 



(i = l~5). (10) 



Here it should be noted that in order to write the above 
expression an extra interaction term is virtually added 
to Hamiltonian (2) for technical convenience as 



Ho(&i,crj,(Jk, <Ji) 

— * Wo(<7j, (Tj, Uk, <Jl), 



<JlOi<Jj 



(11) 



and V3 is, however, put to zero just after differentiation 
with respect to L 3 in eq.(10). 

3. Static Properties 

As shown in I, in order to obtain dynamical suscep- 
tibility x(w) as a linear response to an external field 
E, equilibrium values of order parameters are required 
at each temperature. Since the equilibrium state is 
more easily obtained from the cluster variation method 
(CVM) 9 ) rather than from the stationary solution of the 
time evolution equations(lO), we apply the tetrahedron- 
cactus approximation of the CVM to the present sys- 
tem. 10 > The variational free energy G is obtained by 

G = U-TS, (12) 

where the internal energy U is given by 

(13) 

and the entropy S is given by 



= -V 2 m 2 - V^rrii - V 5 m 5 - 2fi d Em 1 - 2A//.1. 



Nk B i 



= Tr[pi(CTi)lnpi(c7i)] 



- Tr [p 4 , Uj , a k , a t ) In p 4 (a t , Uj , a k , a t )] . 

ijkl 



(14) 



Here c + , c_ , c , d+ , cL and c 2 shown in Fig. 2 are defined 
through p 4 (ai , a ,a k ,ai) as 

c + = ^(l+m 1 +m 2 + m 3 + m 4 + m 5 ), 



c- = 24(1 - mi + m 2 - m 3 + m 4 + m 5 ), 



(15) 



C ° = ^4"(l + m 4-m- 5 ), 

d+ = ^(1 + mi/2 - m 3 /2 - m 4 ), 

d- = - mi/2 + m 3 /2 - m 4 ), 

C2 = -^(1 - m 2 + m 4 + m 5 ), 
with a normalization condition 

c++c_+4co + 4(d + +d_) + 2c 2 = l. (16) 

Since it is sometimes more convenient to use 
c+, c_, c , d+, d_ and C2 instead of mi ~ m 5 especially in 
thermostatic discussions, we rewrite the variational free 
energy in terms of c+, c_, c , d + , d_ and C2 as 

G 

— = 4e c + 4ei(d+ + d_) + 2e 2 c 2 

— 2X tM 2 (p+ -p-f- 2^ d E( P+ - p.) 

— k B T(2p + \np + + 2p_ lnp_ — c + lnc + — c_ lnc_ 

— 4c In c — 4d + In d + — Ad_ In d_ — 2c2 In C2) 
+ s(l - (c+ + c_ + 4c + 4(d+ + d_) + 2c 2 )) , 



with 



Pi(+1) = p+ = c + + 2c + 3d + + d_ + c 2 , 
Pi(-l) = p- = c_ + 2c + d + + 3d_ + c 2 , 



(17) 



(18) 



where s is the Lagrange multiplier to make all the state 
variables c+, c_, co, d+, d_ and C2 independent. With 
A = this expression of the free energy reduces to that 
of Ishibashi. 3 ) Under a uniform electric field the thermal 
equilibrium state is obtained from the minimum condi- 
tion of the free energy: dG/dc + = dG/dc- = dG/dco = 
3G/dd+ = dG/dd- = 8G/dc 2 = as 



co 



Vo 



c+ = (hA) 



4ryo + 2772 + ((hA) 2 + (hA)- 2 ) + A m (hA + (hA)- 1 ) ' 
2 f Co' 



(hA)~ 





c \ 




, C 2 = 7?2 — 




\VoJ 



d+ = (hA)m[-), d_ = (hA)- 1 m (^ 



(19) 
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where 

rji = exp(-/3e 4 ) , (i = 0,l,2) 



A 



1 + m 



cxp(2Dm) . 



(D = /?A M d j (20) 



1 — m 
/i = exp(/3fidE) . 

Further, the polarization m is determined by the relation 

m = p+ -p- = c + - c_ + 2{d + - d-) (21) 

In order to investigate properties of the phase tran- 
sition we make the Landau-type variational free energy 
G(m, E, T) as a function of the electric polarization m 
by substituting (19) into eq.(17): 

G(m,E,T) 
Nk B T 



2Dm 2 - In ■ 



m 



-ln((^ 2 + A" 



4?7o + 4?7i (A + A 1 ) + 2rj 2 ) 

(22) 



It is noteworthy that G(m,E,T) has a property 
dG(m,E,T)/dm = at the thermal equilibrium state. 
Now, we expand G(m, E, T) in terms of order parameter 
m up to the 4-th order: 



G(m,E,T) r 1 „ , mN o 1 „ , mN A 



mE (23) 



4// d (l + 2J>)(l + T ?1 ) 
Tk B T 

where T, A 2 (T) and Aa{T) are given by 
T = 1 + 2770 + 4771 + ry 2 , 

.4 3 (T) = 2(1 + 21?) (l- 2 ™ 1 + '''> ). 
^m-^ + 20 + >5)f4(l + 2D)'(l + m )« 



1 - 2L> 

r 



--(2(l + ?7l ) + (4 + %)(l + 2^) 3 )) 



(24) 



From the view point of Landau's phase transition theory 
the order of the phase transition is classified as follows, 
(i) The phase transition undergoes the second order tran- 
sition at T if A 2 (T ) = and A 4 (T ) > 0. (ii) The phase 
transition undergoes the first order phase transition at 
T c (> T ) if A 2 (T ) = and A 4 {T ) < 0. The boundary 
between the first and the second order is called the tri- 
critical point T t (= T ) when A 2 (T ) = and A 4 (T ) = 0. 
Based on eq.(24) we can get the phase diagram in Fig. 3. 
Since an increase of the phenomenologically introduced 
dipole-dipole interaction strengths the mechanical ferro- 
electric order effect, the ferroelectric transition becomes 
possible before the system comes to an instability point 
when the temperature is decreased. This is the reason 
why the dipolc-dipole interaction induces the first order 
transition instead of the second order transition. As is 



•J* 




tyfd/ 



So 



eJeo 

Fig. 3. The phase diagram of the Slater- Takagi model with dipolc- 
dipole interaction. The solid line stands for the second order 
phase transition line and the dotted one stands for the first order 
phase transition. 



shown in Fig. 3, with ei/en being decreased for each A, the 
transition temperature is decreased because more eas- 
ily thermal excitation of protons favors the paraelectric 
phase. When A 2 (T ) = for A^/e = (Slater-Takagi 
model), A^Tq) — 771/(1 + r)i) is positive definite. That 
is, the Slater-Takagi model always undergoes the second 
order phase transition. Especially, in the Slater model 
(A = 0,£o > 0,£2 > £1 — ► 00) the phase transition takes 
place at the tricritical point. Therefore, it is sometimes 
said that the transition derived by the Slater theory has 
the nature of both the first and the second order phase 
transitions. On the other hand, as is seen in eq.(24) the 
first order phase transition is caused by the finite value 
A/Zd 2 /£o > 0. As A^d 2 /£o is increased, the first order 
transition region becomes larger and the transition tem- 
perature T c is also increased. However, when the dipo- 
lar term becomes dominant, the second order transition 
region again widens because of the phenomenologically 
induced dipolar interaction. 

In order to fix the values of unkown parameters we 
utilize the experimental data on the temperature depen- 
dence of spontaneous polarization 12 ) and choose the pa- 
rameters such as £i/£o = 9.48,£2/£o = 36.92, A^/eo = 
1.095 in the case of first order phase transition. This 
parameter set shows a ratio of spontaneous polariza- 
tion jump to saturated spontaneous polarization about 
0.4. As a comparison with the case of the first order 
phase transition we choose the same parameter set as 
£i/£o = 9.48,£2/£o = 36.92 except A/i 2 /eo = in the 
case of second order phase transition. In the following 
this set of parameters is used in the present calculations. 

Though without dipole-dipole interaction term A = 
we could study the static properties analytically thanks 
to explicit expression of spontaneous polarization, 3 ) with 
A 7^ we cannot obtain analytical spontaneous polarl- 
ization anymore. However, we can deal with these simul- 
taneous equations numerically by making use of the nat- 
ural iteration method (NIM) 13 ) for the CVM. Through 
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Fig. 4. The temperature dependence of spontaneous polarization. 
The solid line (first order) and dotted one (second order) are 
shown. In the first order case the spontaneous polarization jump 
at T = T c is about 0.4. 
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Fig. 6. The uniform susceptibility versus temperature. The solid 
line (first order) and dotted one (second order) are shown. 
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Fig. 5. The entropy versus temperature. The solid line (first or- 
der) and dotted one (second order) are shown. 

the NIM calculation we can obtain the physical quanti- 
ties such as the spontaneous polarization P = 2/ZdWo 
per PO4 and the entropy S for each case of the first and 
the second order transition (Fig. 4 and Fig. 5). These 
numerical results will be used to calculate dynamical sus- 
ceptibility 

The static uniform susceptibility Xstat of the system is 
also obtained as the linear response Am from m induced 
by a uniform field E from the equation of state (21) as 

_ 2Nn d Am _ 2iV/z d 2 1 

X " ' ~ ' " *? (^ + 2^ 

(25) 

with 

A = \j , +m ° exp(2Dm ), 
V 1 - mo 

a = (A 2 + A ~ 2 + 4r7o + 4t7i(A, + AT 1 ) + 2 m y\ 
X = 2a (V + AT 2 + (A + A^- 1 )^) - 2m 2 , 

(26) 



where mo is given from eq.(21) by 

m - a (A 2 - A^ 2 + 2m(A - AT 1 )) . (27) 

When D = 0, this result reduces to our previous result 
for the Slater- Takagi model without dipole-dipole inter- 
action. The numerical result of uniform susceptibility 
Xstat is shown in Fig. 6. 

4. Dynamical Susceptibility 

By making use of the numerical results obtained in the 
previous section, let us calculate dynamical susceptibility 
in this section. When an external electric field E(t) = 
Eexp(iu)t) is applied, we assume that 

m t (t) = m° + x t (u)Eexp(iLut) (i = 1 ~ 5) (28) 

where the required dynamical susceptibility is given by 
x(oj) = fidXii - 1 ) I '4 since it is the response of electric po- 
larization to a weak electric field, and m^'s are equilib- 
rium order parameters in the absence of an external field. 
It should be noted that m®(i = 1 ~ 5) are obtained nu- 
merically by the NIM calculation. By substituting these 
relations (28) into (10) we get the following algebraic 
equations up to the linear ordre to external field: 

(iwTuJ + M) x(w) = b, (29) 

where x(w) = (xi(u>) x 2 (w) X3M Xi{w) XsH)* and 
I is the 5x5 unit matrix. Since the elements of 5 x 5 
matrix M and a 5 x 1 column vector b are a little lengthy 
quantities expressed in terms of only 770,771,772 and D, 
we omit them here. After numerical calculations by the 
Gaussian elimination method for a fixed lutq we get the 
dynamical susceptibility over all the temperature region. 
With the relaxation times T,(i = 1 ~ 5) and the intensity 
coefficients x% the dynamical susceptibility per proton 
x(w) can be written as 

i—l 
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where the five relaxation times n are obtained by a di- 
agonalization of M in terms of a matrix U as 

(UMU- 1 ) t3 = —Sij , (6 i:i : Kroncckcr's delta) (31) 

n 

and the intensity coefficients x* are given by 



nUijbj . 



(32) 



Especially in the paraelectric phase we note that, un- 
der the inversion of an external field E, the order param- 
eters mi and m.3 are changed into —mi and -m3, respec- 
tively, while m-2, m.4 and m.5 are invariant. Accordingly, 
mi and m 3 are long range order parameters responding 
linearly to external field E, while m-2, m.4 and are 
short range order parameters responding quadratically 
to the field. Thus in the paraelectric phase eq. (29) 
is reduced to a closed algebraic equation in mi and 013 
space and as a result the dynamical susceptibility (30) 
can be written down as 



x 



+ ■ 



X 



where 



1 + iujt + 1 + 'mot. 



(7+ = 7-l,T_ = T3) , 



(33) 



p = 2 1 + 



g = 8^ 2 (2V^+V^) 1- 



2^ 2 (1 + 2D) J , 



2(l + 2£»)(l + r?i) 



X 



with 



(t+ - t_)t 

%j3n\K 2 T + T- 
(t+ -t_)t 



r r 

8(l + 7 ? i)(2 v ^p + y^) r_ 

r r 



(1 + 2^/* + ^)^ 



(34) 



(35) 



Next, we plot the real and imaginary part of dynamical 
susceptibility over all the temperature region in the first 
order transition in Fig. 7 and Fig. 8. We see a small 
dip in the real part of x(w)(Fig. 7) consistent with the 
experiments. 4 ^ We notice from Fig. 9 and Fig. 10 that 
one of the five modes overwhelms other four modes in 
relaxation time r» and intensity coefficient x*, though 
any relaxation time r, is finite at T = T c in the first 
order phase transition. That is the reason why the small 
dip structure appears around the transition temperature 
in the real part of the dynamical susceptibility. On the 
other hand, as shown in I, in the Slater- Takagi model 
the phase transition always becomes the second order 
and the most dominant mode shows a critical slowing 
down which leads to the vanishing of the real part of 
dynamical susceptibility at the transition temperature. 
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Fig. 7. The real part of dynamical susceptibility in the case of first 
order phase transition. 
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Fig. 8. The imaginary part of dynamical susceptibility in the case 
of first order phase transition. 
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Fig. 9. The five relaxation times versus temperature in the case of 
first order phase transition. In the paraelectric phase 72, 74, T5 
are omitted because any contribution to x( w ) from them do not 
exist. 
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Fig. 10. The five intensities versus temperature in the case of first 
order phase transition. One of the intensities x 1 becomes dom- 
inant and x 2 ~ X 5 have almost zero value in comparison with 
X 1 over all temperature region. 
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Fig. 11. The frequency dependence of light scattering intensity in 
the case of first order phase transition. 



The imaginary part of dynamical susceptibility is also in 
qualitative agreement with the experiments. 4 ) 

5. Light Scattering Intensity 

Let us calculate the light scattering intensity assuming 
the first order transition with the same previous set of 
parameters and compare it with the experimental results. 
The light scattering intensity I{uo) is expressed in terms 
of the imaginary part of m the high temperature 
approximation as 



I(uj) cx 



Im[xM] 



(36) 



The numerical results of I(u) are shown in Fig. 11. 
When the transition temperature is approached from 
high temperature side above T c , the height of a peak 
increases and the width of the peak decreases around 
lu = 0. 7 ) However, since it is assumed to be the first or- 
der transition, the height of the peak becomes the largest 
at the transition temperature but still finite. Moreover, 
as the temperature is decreased across T c , the height of 
the peak around lu = is decreased drastically. The cen- 
tral mode at u> — is caused by the dipolar fluctuations 



of relaxation modes in the present model. Further, since 
from the imaginary part of cq.(30) it is easily obtained 
as 

/oo poo i 5 

these behaviors are easily understood as follows. Since 
the total integral value of scattering intensity over lu is 
proportional to the sum of intensity of relaxation modes 
at each temperature T, the largest intensity over u> is 
reached at the transition temperature T c but does not 
diverge in the present first order transition. These behav- 
iors agree qualitatively with the experimental data. 11 ) 

6. Conclusions and Discussions 

In order to explain a small dip structure around the 
ferroelectric transition temperature in the real part of 
dynamical susceptibility of KDP crystal observed in the 
experiment, we extend our previous calculation in the 
Slater- Takagi model to the Slater- Takagi model with 
phenomcnologically introduced dipole-dipole interaction. 
We applied the path probability method(PPM) to the 
Slater- Takagi model with dipole-dipole interaction in the 
tetragonal cactus approximation which can take care of 
the ice rule characteristic of KDP crystal. With the 
dipole-dipole interaction being increased, the phase tran- 
sition changes from the second order into the first order 
transition. Further, in the present model the phase tran- 
sition shows the order-disorder transition of proton con- 
figurations. We see the dynamical fluctuations of elec- 
tric polarization modes due to proton configuration in 
the dynamical susceptibility which is expressed in terms 
of relaxation times of the electric polarizations with in- 
tensity of each mode. The Slater- Takagi model always 
induces the second order transition and shows a critical 
slowing down as to relaxation time of the relevant mode. 
Then the system can not follow the external electric field 
with finite frequency at the transition temperature T c 
due to the critical slowing down. This is the reason why 
the real part of the dynamical susceptibility gives zero 
at the critical temperature T c . On the other hand, the 
Slater- Takagi model with the appropriate dipole-dipole 
interaction undergoes the first order transition and the 
relaxation time of each mode becomes finite. Thus the 
system can manage to follow the frequency dependent ex- 
ternal field, though the relaxation time of relevant mode 
is relatively long because of the first order transition close 
to the second order transition of KDP. This fact explains 
why the dip structure is observed around the transition 
temperature in the experiment of the real part of dynam- 
ical susceptibility. We also calculated the light scattering 
intensity utilizing the imaginary part of the present dy- 
namical susceptibility. Since our model deals with relax- 
ation modes in the order-disorder transition, we obtain 
only the central peak structure around to = 0. This peak 
structure increases toward the transition temperature, 
though the height of the peak does not diverges even the 
transition temperature in the present first order transi- 
tion. This behavior is in good agreement with the data 
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of Tominaga et al. ir> 

Finally we should note that all our calculations show 
not quantative but qualitatively good agreements with 
experiments which are ascribed to the Slater- Takagi 
model with phenomenologically introduced dipole-dipolc 
interaction. 
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